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We consider the Higgs boson production at high energy hadron colliders in the framework of 
the fc-r-factorization approach. The attention is focused on the dominant gluon-gluon fusion 
subprocess. We show that fcr-factorization gives a possibility to investigate the associated 
Higgs boson and jets production. 



1 Introduction 

The search for the Higgs boson takes important part at the Fermilab Tevatron experiments and 
will be one of the main fields of study at the CERN LHC collider. In QCD the gluon-gluon 
fusion gg — ► H is the dominant inclusive Higgs production mechanism at LHC conditions. In 
this process, the Higgs production occurs via triangle heavy (top) quark loop. 

It is obvious that the gluon-gluon fusion contribution to the Higgs production at LHC is 
strongly dependend on the gluon density xG(x,fi 2 ) in a proton. Usually gluon density are de- 
scribed by the Dokshitzer-Gribov-Lipatov-Altarelli-Parizi (DGLAP) evolution equation, where 
large logarithmic terms proportional to In// 2 are taken into account. However, at the LHC en- 
ergies typical values of the incident gluon momentum fractions x ~ mn/y/s ~ 0.008 (for Higgs 
boson mass ma = 120 GeV) are small, and another large logarithmic terms proportional to 
lnl/x become important. These contributions can be taken into account using Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) evolution equation. Just as for DGLAP, in this way it is possible to 
factorize an observable into a convolution of process-dependent hard matrix elements with uni- 
versal gluon distributions. But as the virtualities (and transverse momenta) of the propagating 



gluons are no longer ordered, the matrix elements have to be taken off-shell and the convolution 
made also over transverse momentum ky with the unintegrated (fey-dependent) gluon distribu- 
tion J-(x,\£ji). This generalized factorization is called fey-factorization. 

In the collinear factorization, the calculation of Higgs production processes is quite complicated 
even at lowest order because of the heavy quark loops contribution. For example, in Higgs + one 
jet production, triangle and box loops occur, and in Higgs + two jet production the pentagon 
loops occur (see^ and references therein). However, the calculations of the Higgs production 
rates can be simplified in the limit of large top quark mass mj — ► In this approximation the 
coupling of the gluons to the Higgs via top-quark loop can be replaced by an effective coupling. 
Thus it reduces the number of loops in a given diagram by one. The large m< approximation 
is valid to an accuracy of ~ 5% in the intermediate Higgs mass range m# < 2mt, as long as 
transverse momenta of the Higgs or final jets are smaller than of the top quark mass (pr <mt)- 

A particularly interesting quantity is the transverse momentum distribution of the produced 
Higgs boson. It is well-known that the fixed-order perturbative QCD is applicable when the 
Higgs transverse momentum is comparable to the m#. Hovewer, the main part of the events 
is expected in the small-py region (pr *C ttih), where the coefficients of the perturbative series 
in a s are enhanced by powers of large logarithmic terms proportional to mm^/p^. Therefore 
reliable predictions at small px can only be obtained if these terms will be resummed to all 
orders. Recently it was shown^ that in the framework of fey-factorization approach the soft 
gluon resummation formulas are the result of the approximate treatment of the solutions of the 
CCFM evolution equation (in the 6-representation) . 

There are several additional motivations for our study of the Higgs production in the fey- 
factorization approach. First of all, in the standard collinear approach, when the transverse 
momentum of the initial gluons is neglected, the transerse momentum of the final Higgs boson 
in gg — ► H subprocess is zero. Therefore it is necessary to include an initial-state QCD radiation 
to generate the py distributions. In the fey-factorization approach the underlying partonic 
subprocess is gg — » H and the fey-factorization naturally includes a large part of the high-order 
perturbative QCD corrections'^. Since the upper gluon in the parton ladder is not included 
in the hard interaction, its transverse momentum is now determined by the properties of the 
evolution equation only. It means that in the fey-factorization approach the study of transverse 
momenta distributions in the Higgs production via gluon-gluon fusion will be direct probe of 
the unintegrated gluon distributions in a proton. In this case the transverse momentum of 
the produced Higgs should be equal to the sum of the of the initial gluons. Therefore future 
experimental studies at LHC can be used as further test of the non-collinear parton evolution. 

In the previous studies BEEI ^ e fe r -factorization formalism was applied to c alcu late trans- 
verse momentum distribution of the inclusive Higgs production. The calculations^^ were done 
using the on-mass shell (independent from the gluon fey) matrix element of the gg — > H sub- 
process and rather the similar results have been obtained. In Ref.^ in the framework of MC 
generator CASCADE^! the off- mass-shell matrix element obtained by F. Hautmann^ has been 
used with full CCFM evolution. In our paper^we have investigated Higgs production at hadron 
colliders using the full CCFM-evolved unintegrated gluon densities. Here in order to illustrate 
the fact that in the fey-factorization approach the main features of collinear higher-order pQCD 
corrections are taken into account effectively, we give theoretical predictions for the Higgs + one 
jet and Higgs + two jet production processes. 

2 Basic formulas 

The effective Lagrangian for the Higgs boson coupling to gluons^ is 




(1) 



where Gf is the Fermi coupling constant, Cr" is the gluon field strength tensor and H is the 
Higgs field. The triangle vertex T^ u {k\, k 2 ) for two off-shell gluons having four-momenta k\ and 
k 2 and color indexes a and b respectively, can be obtained easily from the Lagrangian (1): 

T^(h,k 2 ) = i8 ab ^ (G F V2) 1/2 \k»k\ - (h ■ k 2 ) 9 n . (2) 

The differential cross section of inclusive Higgs production pp — > H + X in the k F — factorization 
approach has been calculated in ^and can be written as: 



da{pp^H + X) _ f a 2 s {ji 2 ) G F y/2 r 2 2 i2 

m H + PT 



dyn J 2887T xix 2 m H s 



cos 2 (A</?)x 

(3) 



x^(x X) kf T , /i 2 M(x 2 , k| r , /z 2 )a?k 2 T dki r ^^, 

where ^4(x, k 2 ^, /i 2 ) is the unintegrated gluon distribution, A(p the azimutal angle between the 
initial gluon momenta kiy and k2y, and the transverse momentum of the produced Higgs bo- 
son is pt = kiT + k 2 T- It should be noted that this process is particulary interesting in 
k F — factorization, as the transverse momenta of the gluons are in the same order as their longi- 
tudinal momenta (~ (10 GeV))^. 



3 Numerical results 
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Figure 1: Differential cross section da/dpT for Higgs + one jet production at y/s = 14 TeV: solid, dashed, dash- 
dotted and dotted curves correspond to win = 125,100,150,200 GeV, respectively (left column). The jet-jet 
azimutal angle distribution in the Higgs + two jet production: solid and dashed curves correspond to the J2003 
set 1 and J2003 set 2 unintegrated gluon distributions, respectively (right column). 

Our prediction for the total cross secrion and the transverse momentum and rapidity distri- 
butions of the inclusive Higgs production at the LHC (yfs = 14 TeV) were done in RefM The 
calculations were fulfiled for four choices of the Higgs boson mass {run = 100, 125, 150, 200 GeV) 
under interest in the Standard Model with default scale fx 2 = m 2 H . We used LO formula for the 
strong coupling constant a s (fj?) with nj = 4 active quark flavours and Aqcd = 200 MeV, such 
that a s (M§) = 0.1232. 

Here we illustrate how fcT-factorization approach can be used to calculate the semi-inclusive 
Higgs production rates. We choose the one carrying the largest transverse momentum, and 
then compute Higgs with an associated jet cross sections at the LHC energy. We have applied 
the usual cut on the final jet transverse momentum | Pjet T 1 1 > 20 GeV. Our predictions for the 
transverse momentum distribution of the Higgs + one jet production are shown in Fig. 1 (left 
column). One can see the shift of the peak position in the pt distributions in comparison with 



inclusive production, which is direct consequence of the | Pjet T" I > 20 GeV cut. To demonstrate 
the possibilities of the fey-factorization approach, we calculate azimuthal angle distributions 
between the two final jet transverse momenta in the Higgs + two jet production process, where 
the kinematical cut |pjetr| > 20 GeV was applied for both final jets. Studing of these quantities 
are important to clean separation of weak-boson fusion and gluon-gluon fusion contributions. 
Our results are shown in Fig. 1 (right column), where a dip at 90 degrees is seen. This dip 
at A<p ~ 7r/2 comes from the cos(Aip) in eq. (3). In the approach presented here, the kx 
of the initial gluons is approximately compensated by the transverse momenta of the jets ^4 
kr — Pjctr and, consequently, Acp ~ Aip. This dip is characteristic for the loop-induced Higgs 
coupling to gluons in the framework of fixed-order pQCD calculations HI Thus, we illustrate 
that the features usually interpreted as NNLO effects are reproduced in the k^— factorizatiuon 
with LO matrix elements. 

However, we see a verylarge difference between the predictions based on the J2003 gluon 
densities set 1 and set 2™, showing the sensitivity to the shape of the unintegrated gluon 
distribution. 



4 Conclusions 



The predictions in the kx~ factorization approach are very close to NNLO pQCD results for 
the inclusive Higgs production at the LHC, since the main part of high-order collinear pQCD 
corrections is already included in the kx~ factorization. In the kx~ factorization approach the 
calculation of the associated Higgs+ jets production is much simpler than in the collinear fac- 
torization approach. However, the large scale dependence of our calculations (of the order of 20 
- 50%) probably indicates the sensitivity to the unintegrated gluon distributions. 
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